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ABS TRACT 


An investigation to explore the possibility of whirl flutter and to find the effect of pitch-flap 
coupling (63) on teetering motion of the DOE/NASA Mod-2 wind turbine is presented. The equations of 
^isjtion are derived for an idealized five-degree-of- freedom mathematical irodel of a horizontal-axis 
wind turbine with a tw o-bladed teetering rotor. The nnodel accounts for the out-of-plane bending 
flM3tion of each bl ade, the teetering motion of the rotor, and both the pitching and yawing motions of 
the rotor s^^rt. Results show that ^e”ltod-2 design is free from whirl flutter. Selected results 
are pra^ nt¥d~ rndicatir^ the effect o f var iations in rotor support dafi^ing, "rotor support stiffness, 
pitching, yawing, teetering, and blade bending motions. 


INTRCDUCTION 

Recent hori zont al-axis wind turtine (HAWT) 
designs such as the DOE/NASA hfc^d-2 wind turbine 
(ref. 1) include flexible towers in order to 
achieve significant weight and cost reductions. 
Experience with prop-rotors has shown that rotors 
with flexTbTe" supports have a potential 
^roeIastrc~instability known as whirl flutter. 
This form of "Instability the Interaction 

of elastic,” darrping, gyroscopic, arxl aerodynamic 
forces. The whirl flutter problem is discussed 
in references 2-7 amor^ others. In whirl 
instaoility, the rotor will precess In a whirl 
TOde with an eve r-iryreasi ng amplitude wh en the 
critical wind speed has been reached. That is, a 
point on the rotor hub will trace a divergent 
spiral as illustrated in Figure 1. The direction 
Qf_ths ^iral rotation can be either the same as, 
or counter to, the rotor rotation. These two 
modes are referred to as fo rward and ba ckward 
whirl HKides, respectively. Contimed operation 
of a wind turbine in the whirl flutter 



Figure 1. - Wind Turbine Rotor in a Forward Whirl 
^tode. 


mode will quickly lead to failure of~th¥ 
supportive structure. This v^irl instability is 
possible regardless of the pfesehce of rotor 
teetering motions or blade out-of-plane bending 
motions. When these n(K)tions are In clu ded, they 
couple wiitt the nations of the supportive 
structure. Then a v^irl instability can occur in 
the whirl modes of the supportive structure 
afxj_/or the rotor. ' 


Most of the current large HAWT systems have 
rotors with two blades. The analysis of wind 
turbines with two-bladed rotors differs 
signific¥^tl7 f^TO*^ that with axisyrmnetric 

rotors. T^ properties of a two-bla ded wind 

"turbine change significantly as the blades” rotate 
from a horizontal to a vertical position. As a 
result, the equations of motion of a two-bladed 
wind turbine system contain significant periodic 
coefficients. 

In order to reduce blade bending loads, teetered 
rotors with pitch-flap coupling (<53) have been 
used in sorr^ HAWT systems. The pitch- flap 
coupling mechanically changes the pitch of the 
blades as the rotor “teeters and thus is 
equivalent to an aerodynamic spring that 
restrains teeterir^ TOtion. The effect of 63 
on rotor motion stability was studied in refs. 6 
and 8. A whirl flutter analysis for a prop -rotor 
on a flexibly runted pylon was developed in ref. 
6. That analysis may also be suitable for 
investigation of whirl flutter in HAWT systems. 
However, since most HAWT systems use rotor blades 
that are long and relatively flexible, the blade 
flexibility ought to be included in the 
formulation of a HAWT whirl flutter analysis. 

Other analyses are available for the 
investigation of whirl flutter in HAWT systems. 
One is the ^C}STAS corrputer code (refs, 9 and 
10). However, it is very complex and uses a 
large amount of computer time. Hence, it is not 
well suited to parametric investigations. 
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Therefore, a simple model, enconpassing only the 
pertinent degrees of freedom, is desired to study 
the possibility of whirl flutter in a flexibly 
mounted HAWT. 

The primary purpose of this paper is to present 
the development of a simple model for exploring 
the possibility of whirl flutter in the DOE/NASA 
Mod-2 HAWT. Secondary purposes are to study the 
effects of pitch-flap coupling, rotor support 
stiffness, and rotor support darrping on the 
response of the Mod-2. 

A five-deg re e-of- freedom mathematical model is 
developed in the Appendix for a flexibly-mounted 
two-bladed teetering rotor. The degrees of 
freedom include the first out-of-plane bending 
mode for each blade, the rotor teetering motion, 
the rotor support pitching motion, and the rotor 
support yawing motion. The developed equations 
that have periodic coefficients are numerically 
integrated in the time domain using a standard 
Runge-Kutta method. 

ANALYSIS NtTHOD 
Mathematical Model 

The mathematical model of a HAWT with two-bladed 
teetering rotor is shown in Figure 2, The rotor 
support is modeled by a rigid pylon of length h 
that is restrained at one end by two sets of 
rotational springs and dampers. These springs 
and dampers represent tower stiffnesses and 
dampings. The restraints allow only pitching and 
yawing motions, ^y» pylon. The 

teetering motion, Y, of the rotor hub with 
respect to the rotating shaft of the pylon is 
also restrained by a rotational spring and damper 
set. The angular velocity, of the rotor is 
assumed to be constant. The out-of-plane blade 
bending deflections are represented by and 
W2. These deflections are, in turn, expressed 
in terms of the normal bending modes and the 
generalized coordinates. Since the blades are 
relatively stiff, only one mode is considered. 
This type of representation of the blade motion 
is referred to as a Rayleigh-type of analysis. 

As a consequence of this approximation of the 
blade motion, there are three degrees of freedom 
for the rotor, one for each blade, and one for 

y 



Teetering HAWT. 


teetering. Thus, with the pitching and yawing 
motions of the pylon, the wind turbine firadel has 
a total of five degrees of freedom. Only the 
out-of-plane bending motion of the blades is 
considered because it couples with the rotor 
teetering motion. Consideration of other notions 
such as tower translation, blade in-plane 
bending, and blade torsion are rot difficult, but 
their inclusion would increase the complexity of 
the analysis. Furthermore, it is believed that 
these other motions do not have much effect on 
vhirl flutter. 

The aerodynamic forces are obtained from strip 
theory based on a quasi-steady approximation of 
two-dimensional, inconpressible, thin airfoil 
theory. The blade geometric pitch angle, which 
consists of the blade built-in twist (pretwist), 
the pitch angle due to pitch-flap coupling, the 
collective pitch angle, and the cyclic pitch 
angle are included in the formulation. Classical 
blade element momentum theory is used to 
calculate the steady induced velocity. 

Coordinate Systems 

Several orthogonal coordinate systems are used in 
the derivation of the equations of motion. Those 
that are common to both the dynamic and 
aerodynamic aspects of the HAWT are described in 
this section. 

1. Inertial system XYZ — The Y-axis of this 
system, shown in Figure 2, coincides with the 
vertical axis of the HAWT tower and is positive 
upward. The Z-axis coincides with rotor axis and 
is positive into the wind. 

2. Hub system X3Y3Z3 — This system is 
fixed to the hub center but does not rotate with 
the rotor. It is parallel to the XYZ system when 
the pod rotations are zero. 

3. Rotor system X4Y4Z4 — This axis 
system is obtained by rotating the hub system 
about the Z3 axis by the rotor position angle ^ 
(=Qt) as shown in Figure 2. 

4. Blade system xpytjZp — This axis 
system is obtained by rotating the rotor system 
about the X4-axis by the rotor teetering 
angle y. The yp axis is aligned along the 
blade quarter chord points and is also assumed to 
be the blade elastic axis. The xp and zp 

axes are also shown in Figure 3 along with the 
various blade element angles, relative 
velocities, and resultant aerodynamic forces. 



Figure 3. - Blade Element Velocity and Force 
Vectors for a wind Turbine Rotor, 
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Computer Code 


The equations of motion developed in the Appendix 
have timewise periodic coefficients. The 
stability of a HAWT rmjst be determined by 
numerically integratir^ these equations or by 
-using Floquet-Liapunov theory. To this end, a 
compKjter program called ASTER5 (Aeroelastic 
Stability of a TEetering ^otor wTth 5 degrees of 
was" written to numerically Tntegrate 
these equations. The ASTER 5 program was first 
verified by several special cases obtained from 
ref. 6. the program was then used to investigate 
the possibility of whirl flutter in the DOE/NASA 
Mod-2 HAWT and the effect of variations in some 
of the Mod-2 parameters on its response. 


The ASTERS conpu ter program was ^written in 
FORTRAN Iv, The input includes the radial 
distributions of blade chord, twist angle, mass, 
3hd first cut-of-plane bending mode; equivalent 
-inertia, stiffness, and dampir^ constants for the 
pylon; and aerodynamic data. The input allows 
■part-span pitchable blades with pitch- flap 
t otjpl lng and cyclic pitch. The program uses a 
st and ard sttorout ine call ed DV ERK, which solves a 
-systam ~of f irst^x^er~dif equations with 

a Runge-Kutta method based on Verners fifth- and 
sixth-order pair of formulas. 


RESULTS AND DISCUSSIOl 

To verify the ability of the ASTER5 program to 
correctly predict whirl flutter, several cases of 
a prop-rotor, which was analyz^ in ref, 6, were 
evaluated. The parameters for the prop-rotor are 
presented in Table I. Results are presented for 
two typical cases. In one the prop-rotor 

response exhibits whirl flutter, while in the 
other case it is stable. The whirl flutter case 
in which the pitching and yawing frequencies are 
2.3 Hz and 5.0 Hz respectively, is shown in 
Figure 4. When the pitching frequency is raised 
to 3.3 Hz by increasing the pitch spring 

stiffness, the prop-rotor becorr^s stable, as 

shown In Figure 5. For comparison, the envelopes 
of the pitch motion amplitudes for the 
corresponding cases calculated in ref. 6 are also 
Indicated in Figures 4 and 5. It is evident from 



Figure 4. - Response of Prop-Rotor in 
Whirl Flutter ^de. 


these figures that there is agreement between the 
results of ref. 6 and the ASTER5 program. Thus, 
the ASTERS program is capable of predicting whirl 
flutter. The quantitative differences evident in 
U^ese figures may be due to differences in 
airfoil data and/or initial conditions. It 
should also be noted that ref. 6 does not account 
for the blade out-of-plane bending motions as 
does ASTERS. However, the blade frequency is 
assumed to be high for the input to ASTERS, and 
thus has a negligible effect on stability. The 
steady state pitch deflection, evident in 
Figure S, is due to the gravitational moment of 
the rotor, which is added to the pitching moment 
only for these verification cases. 



TIHf, t. 


Figure 5. - Response of Prop-Rotor in 
Stable Mode. 


The DOE/NASA Mod-2 HAWT was modeled to 
investigate the possibility of whirl flutter. 

The parameters for Mod-2 are presented in Table 
II, The response of the Mod-2 was calculated 
with the ASTERS program. A baseline reference 
xase of the Mod-2 parameters without structural 
damping was considej^ed for an initial evaluation 
of its stability. The~ results of this case, 
given in Figure 6, show that the pitch, yaw. 
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teeter, and blade cyclic bending motions are 
neutrally stable. However, when a small amount 
of structural damping for the pitch and yaw 
motions is Included, all motions are damped out 
as shown in Figure 7, Since damping exceeding 
this amount is expected in the actual system, it 
is concluded that the baseline Mod-2 is free from 
whirl flutter. 
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Figure 7, - Response of Mod-2 Baseline Case with 
Structural Damping 


To study the effect of pitch- flap coupling on the 
response of the baseline Mod-2, several cases 
were calculated for values of % from -40° to 
+40°. The results indicate that only the blade 
cyclic bending motion as measured by Oq is 
affected by variations in Figure 8 shows 
the change in maximum amplitude of Qq with 
The results indicate that positive ^ has an 
adverse effect on blade cyclic bending motion. 



PITCH OAP COUPLING. deg. 

Figure 8. - Effect of Pitch-Flap Coupling, ^3, 
on Blade Cyclic Out-of-Plane 
Bending Motion, 

Other parametric studies were made to explore the 
possibility of whirl flutter over wide ranges of 
pylon spring stiffnesses, pylon dampings, rotor 
rotational speeds, and wind speeds. Some 
selected results of these studies are presented 
in Fibres 9-12. The possibility of whirl 
flutter can exist for Mod-2 if the yaw or pitch 
stiffness of the pylon were substantially 
reduced. For example, Figure 9 shows the 
response of Mod-2 when the yaw stiffness is 
decreased to 6.6% of its baseline value while the 

other parameters remain the same. These results 


•indicate vrfiirl flutter by the unstable response 
of the yaw and teeter motions. When the pylon 
pitch stiffness is also reduced to 7.3% of its 
baseline value such that the pitch and yaw 
frequencies are equal (^(j)x=^y=3.665 Hz), 
then the response of the pitcn motion is also 
unstable as shown in Figure 10. The whirl motion 
of the pylon for this case is best illustrated by 
a cross-plot of the pitch and yaw motion in 
Figure 11. The figure shows that the system is 
in a forward whirl mode. From these results, it 
can be concluded that the stability of a HAWT is 
hi^ly dependent on the rotor support stiffnesses. 


As demonstrated earlier, the stability of a HAWT 
is sensitive to the presence of structural 
damping. To further illustrate this fact, a 
nominal amount of damping (^4>x=^^y=.04) was 
added to the unstable case of Figures 10 and 11. 
The results, shown in Figure 12, indicate that a 
reasonable arrcjunt of structural damping has 
stabilized all motions of a previously unstable 
system, 
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Figure 9. - Response of Mod-2 with Reduced Yaw 
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Figure 10. - Response of Mod-2 with Reduced 
Pitch and Yaw Stiffness 
(‘^x=“4’y=3.665 Hz). 
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3. Positive 63 has an adverse effect on 
cyclic blade out-of-plane bending nations for the 
Mod-2 design, whereas negative 63 has little 
effect. 

4, Reduction in rotor support stiffness or 
structural danping increases the possibility of 
w^iirl flutter. 
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CONCLUSIONS 

An investig ation was conducted to explore the 
possibility of whirl flutter in a large HAWT. A 
five degree-qf- freedom mathematical model and its 
associated' computer program were developed and 
ve^fied. The program was used to study the 
possibility of whirl flutter in the DC£/NASA 
Hod"-? wind turbine and the effect of parametric 
variations in pitch-flap couplir^, rotor support 
stiffnesses, and structural damping on its 
respcnsa. _ Based on these limited studies, the 
following conclusions were obtained. 

1. The ASTER5 program is capable of 
predicting whirl flutter for two-bladed teetering 
rotor systems, 

“' 2 . baseline design of the Mod-2 HAWT is 
free of whirl flutter. 
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^MENCLATURE 


c blade chord length 

Cp, Cw, damping coefficients of rotor 

teetering, blade out-of- 
^ plane bending, pylon pitch and yaw 

motions, respectively 

D profile drag per unit length of 


blade element 
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h 

H 

Hi, H2 

Ib 

^Px’ ^Py 
Kb 

Kq f Ky^ 

Nix’ 

L 

Mb, 

Qc » 

^W]^> ^W 2 

r 

R 

Sbi 


pylon length 

X 4 , Y4, Z4 

rotor coordinate system 

total rotor shear force, Eq. (A20) 

a 

blade angle of attack 

rotor shear force per unit length 

Y 

rotor teeter deflection angle 

of blade 1 and 2 , respectively 


pitch-flap coupling angle Eq. (/ 

mass moment of Inertia of the blade 

^w» 

critical damping ratios, 

defined in Eq. (A9) 


Eq. (A13) 

mass moments of inertia of 

9 , 9t, 00 

blade pitch, twist, collective 

the pylon about the X and Y axes 

^IC, ^15 

pitch, and cyclic pitch angles 

half of the teeter spring stiffness 

P 

air density 

effective blade spring stiffnesses 
defined in Eq. (All) 


aerodynamic inflow angle 

pylon spring stiffnesses 

^x» ^y 

pylon rotational deflections 

circulatory lift per unit length of 


normalized blade mode shape 

blade element 


rotor position angle 

mass properties of the blade 
defined in Eq. (A9) 

wb 

blade natural frequency 


^w 

blade out-of-plane bending 

cyclic and symmetric coordinates 
for blade out-of-plane bending 


frequency 

motions defined in Eq. (A9) 

^^x 

pylon frequency ) 

generalized coordinates 

for out-of-plane bending motions of 

a)(J)y 

pylon frequency (=M^y/I,py ) 

blades 1 and 2 

U 

rotor rotational speed 

radial distance along blade elastic 
axis 

(*) 

time derivative 

radial length of blade 

{ } 

column matrix 

mass property of blade defined in 

[] 

square matrix 

Eq. (A9) 

TABLE I. 1 

PARAMETERS FOR PROP-ROTOR OF REF. 6 


t 


time 


T total rotor thrust force, Eq. (A20) 

also kinetic energy 

Tl, l 2 rotor thrust force per unit length 

of blade 1 and 2 


U 

Up, Uf 


V 

V 


w^, W 2 


resultant aerodynamic velocity, 
also potential energy 

coirponents of U, Figure 3 

induced velocity 

wind velocity 

out-of-plane bending deflections of 
blade 1 and 2 


^b> Ybj ^b blade coordinate system 
X, Y, Z inertial coordinate system 

X 3 , Y 3 , Z 3 hub coordinate system 


Air velocity, V 77.1 m/s 

Rotor 

Radial length, R 3.505 m 

Rotational speed, 320 RPM 

Pitch-flap coupling, 63 20° 

Teeter spring stiffness, 0 

Teeter motion damping, 2 0 

Blade 

Mass properties 




32 

Sb? 

Ib 

Stiffness, (KQ+Ky,) 
Damping, 

Airfoil 
Chord, c 

Twist distribution, 
0 <r/R ^ .45 
45<r/R <1.0 
Collective pitch, 
Pylon Properties 

Inertias, Ip = Ip 
Stiffness, ^ ^ 

K<j)x 

K(|)y 

Damping, 

Length, h 


26.09 kg 
16.96 kg-m2 
42.17 kg-m2 

118.9 kg-m2 
.3821x10^ N-m/rad 

0 

NACA 0015 
.2794 m 

9t(r) 

.677-1. 21 7r/R rad 
.419(.75-r/R) rad 
9q .74 rad 


21,60 kg-m2 

29.71x10^ N-m/rad 
140.3x10^ N-m/rad 
.04 

1.143 m 
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TABLE II* PARAtCTERS TOR DOE/NP 

Wind v/elocity, V 
Rotor 

Radial length, R 
Rotatinal speed, 

Pitch-flap coupling, 63 
Teeter spring stiffness, 2Kt) 
teeter nration darnping, 2 cp 
Blade 

Mass properties 


32 - 

3.368x10^ 

Stiffness, (Kq+K^) 78.82x10^ 

Damping , Cw 0 

Airfoil NACA 23018 

Chord distribution, c(r) 

. 1542 ^r/R^ *3455 3.319-8.429( ,3455-r/R) m 

.3455 <r/R^ 1.0 1. 436+2. 877( 1-r/R) m 

Twist distribution, 

. 1542 < r/R < . 27 . 03459- . 155( . 27-r/R ) rad 

,27<r/RS1.0 -.0698 + .143(l-r/R) rad 


Collective pitch, ^ 
- .1542< r/R ^ .7006 
;7TO6sr/R^1.0 

Pylon 

_ Inertia, 

^Px 


Ipy 

tifrness, 


«>y 

LergU-., h 


0 

-.05236 rad 


6.115x106 N-m2 
.6210x106 N-rt|2 

6.183x1q 9_ N-fti/rad 
3.140 x1o 9 N-m/rad 

0 

0 

7.3152 m 


HAWT 

where 




m/s 


r 

1 

0 

0 

m 

RPM 

II 

0 

cos d) 

-sin 



0 

L. 

sin * 

X 

cos 0 

^<g 

kg-m2 

kg-m2 

N-m/rad 


r 

cos 0 0 

y 

0 1 

-sin 0 

L y 

sin 

0 

cos 4^ 


(A3) 


The position vector of a point on blade 1 can be 
written in the X 3 Y 3 Z 3 axis system as 


r, = [xJ[T ] { 


~\ 


wi(r,t)j 


(A4) 


»*iere 


M- 


^OS 

sin 

^ 0 


-sin ip 
cos ^ 

0 


0 

0 

ly 


APPENDIX 


0 0 




(A5) 


Derivation of Equation of Motion 


Tne mathematical nradel of a horizontal axis wind 
turbine is shown in figure 2. The degrees of 
freedom and the required coordinate systems are 
described in the main body of this paper. The 
^uations of nration, herein, have been derived by 
usif^ the Lagrangian approach. This formulation 
requires expressions for the position vectors of 
arbitrary TXDints on the pylon and the blades. 

These expressions are obtained with the aid of a 
series of rotations. The order of the rotations, 
illustrated in Figure 2, is Y. 

The position vector of a point on the pylon axis is 


r 

P 


' 0 ^ 




(Al) 


K] = 


^0 


COS Y -sin Y 
sin Y cos Y . 


Combining Equations (A2) and (A4), the position 
vector of a point on the axis of blade 1 expressed 
in the XYZ axis system is 


r - 


where wx(r,t) Is represented by a single 
elastic blade mode and is 






(r.t) 


(A 6 ) 


and that of the hub-pylon axis junction point is 

Co' 

Ch ' P'<}>xKT4,y] 


(A2) 


w^(r,t) = ^(r)q^^(t) (A7) 

The position vector of a point on the axis of 
blade 2 is obtained from Equations (A 6 ) by 
replacing wi, Y, and 4^ by W 2 , -Y, and 4^+1i 
respectively. 


207 



where 


The total kinetic energy of the pylon and the 
rotor is formed from the position vectors given by 
Equations (Al) and (A6) and is given by 

T - y 

+ + Ijj(l - COS 2}p)] 

+ - 21^y^ + 


‘'w' 

\ “ jf" ^c'^’^(r)dr 


(All) 


+ + 2Mb^(q; + q2) + sin 2^ 

+ + ^b’'^ '1^ 

+ ^4’x’I^Ib \l/ + 4$^q^Sj,^ cos Ip 

- 40^0 jij^i).^ + S^^q^ cos Ip + I^Y cos ipj 


(A8) 


+ 4S, 


>X^^c} 


where 




S ■ ^Py 


dr 


dr 


r dr 


(A9) 


’^•£- 

h ■ v' 

S 

\ '£ v“ '■' 


■*■ 


**" 

The quantities ^p>^ and ^py are the pylon 
inertias about the X and Y axes, respectively. 

The potential energy of the pylon and the rotor 
can be written as 


« - iK-^x + n/y * 


+ 2(K^ + K^)(q3 + q: 


The dissipation potential for the pylon and the 
rotor can be written as 


“d ■ + 2C^Y^ 


+ 2C 


where 


-(^C ^s3 


(A12) 


~ ^^lpy^lpy“lpy 


2^^b“Y 




(A13) 


By substituting Equations (A8), (AlO), and (A12) 
into Lagrangian equations of the form 


J- fjl.\ _ JL + JS. j. ^ . n 

yhi) 3q^ 3qj^ 3q^ '^1 


(A14) 


the following equations of motion for the wind 
turbine model are obtained 


where 


k} 


q} + 

[C]{q} + [K]jq| - |Q} 




**c 


Qwi “ Qw 2 

V 


Hy y 



M<P^ 





(A15) 


[I] 


0 

0 

0 

0 


0 0 0 0 
2%2 ^ * 

2Sbj 21b ^ 

2Sbi coi it 2Ib coi + IbU + cos 2i|J) Ib sin 2’)' 

2Sbi sin 2Ib •In ij> Ib sin 2i|/ Ifa^^ ‘ 

(A16) 


W 


(AlO) 


2C 0 0 0 • 0 

w 

0 2C^ 0 •!“ 'I' *®bi^ ^ 

0 0 2C^ -4lb^ •!“ 'f* ^Ib^ ^ 

0 0 0 - 21^^ sin 2(Jr 2Ib^(l + cos 2 ^>) 

0 0 0 -21bfi(l - cos 2ijii) sin 


w- 


0 

0 

0 


> 0 0 0 0 
2CK* + K^) 2Sbjfl^ 0 0 

ZSbjfl^ 2Kb + 2Ib«^ 0 0 

2Sbjf^^cos Zlb^^cos K^ 0 

2Sbjft^slti 2lbR^«in fp 0 K^y 
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where, from Figure 3 and Equation (A 6 ), the 
following expressions are obtained 


The r^xt step is to obtain expressions for 
Qwx, Q«r2, My , H(>x, and l^y. The 
expressions are derived from the virtual work of 
the aerodynamic forces, which can be written as 




■ij' 


(Fa ’ 6 r)dr 


(A17) 


where is the aerodynamic forc£ vector per 
unit ler^th. The components of are 
illustrated in Figure 3 from ¥^ich one can write 


K\ 


blade 1 


(A18) 


By using a similar expression for blade 2, 
_sUjstitutir^ it with Equations (A 6 ) and (A18) into 
Eq~uatxo n (A 17), and neglecting several higher 
order terms believed to be unimportant, one obtains 




9-X-. 

/ T,$ dr 



h- 

/ dr 

Ad ^ 




' (T| - T, 

) 


,)r dr 


(A19) 


-Hh sin Ilf + Thy cos + cos tp 



Hh cos ip + Thy 


sin Jp + sin \p 


where 



(Hj - H 2 )dr 
(Tj + T^)dr 


(A20) 


u - 



+ r<^y sin ^ + ry + + v - - 


“ V cos sin cos - V sin <p^ sin 
+ sin 4 / - h4»y cos 4^ + rQ cos y (A22) 

a * - 6 

-1 tv 

ip - tan(— ) 

6 - Sq + 6 ^(r) - y tan 6 ^ 

+ 0, cos li* + 9, sin 4^ 

In the derivation of above expressions for Up 
and Uj, several higher order terms, believed to 
be unimportant, are neglected. Also, circulatory 
lift, produced by the angular velocity of the 
local blade section about the yj;j-axis due to 
blade out-of-plane bending and rotor teetering, is 
neglected in the expression of Equation (A21). 

The values of C[_ and Cq are nonlinear 
functions of x and are calculated from airfoil 
data. 

By resolving L and D in Figure 3, the expressions 
for Ti and Hi are 

Tj = -L cos 4> - D sin (p 

_ .. _ . (A23) 

» -L sin (|) + D cos <P 

The expressions for I 2 and H 2 are the same but 
the values of L, D, and are obtained by 
replacing y, 4 ^, and Qwi by -y, 4^+ir, and 
^W 2 , respectively in the expressions of 
Equation (A22), 


“The expressions for circulatory lift and profile The induced velocity v in Equation (A22) is based 
drag per unit length can be written as on classical momentum theory and is 


L - pU^cCj^(a) 

D - pU^cCjj(a) 


(A21) 


V * 


-V + 



-t- 2 T/pttR^ 

2 


(A24) 
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QUESTIONS AND ANSWERS 


D.C. Janetzke 


From: A. Wright 

Q: What type of failures would result from severe whirl flutter? 

A: Fatigue failure or ultimate limit load failure. 

From: Bill Wentz 

Q: How do you increase structural damping in design? 

A: I don ^ t know. 

From: J.A. Kentfield 

Q: What magnitude of structural damping can be expected in the pylon of MOD-2 or 

similar machines? 

A: The damping applied to the pylon in the model represents the equivalent damping of 

the entire rotor support system which includes the pod and the tower. The Mod-2 
welded tower damping is about 2% of critical damping. 

From: Mr. Doman 

Q: What influence has the absence of tower bending modes on results? 

A; The tower bending modes are represented by the pylon support stiffness , 

From; P. Anderson 

Q: What time step size was used in the integration process? Have any sensitivity tests 

been carried out to optimize step size? 

A: 1, An initial time step equal to ZS steps per rotor revolution was arbitrarily ^ 

chosen. The integration process could change the step size within the initial 
size as needed for convergence, 

2, Several other step sizes were used^ but no attempt was made to optimize the 
size. 
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